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ABSTRACT 

We present an analysis of the X-ray spectral properties observed from 
black hole . candidate (BHC) binary SS 433. We have analyzed Rossi X-ray 
Time Explorer ( RXTEJ) data from this source, coordinated with Green Bank 
Interferometer /RATAN-600. We show that SS 433 undergoes a X-ray spectral 
transition from the low hard state (LHS) to the intermediate state (IS). We show 
that the X-ray broad-band energy spectra during all spectral states are well fit by 
a sum of so called “Bulk Motion Comptonization (BMC) component” and by two 
(broad and narrow) Gaussians for the continuum and line emissions respectively. 

In addition to these spectral model components we also find a strong feature 
that we identify as a ”blackbody-like (BB)” component which color temperature 
is in the range of 4-5 keV in 24 IS spectra during the radio outburst decay in 
SS 433. Our observational results on the ’’high temperature BB” bump leads 
us to suggest the presence of gravitationally redshifted annihilation line emis- 
sion in this source. In fact this spectral feature has been recently reproduced in 
Monte Carlo simulations by Laurent and Titarchuk. We have also established 
the photon index saturation at about 2.3 in index vs mass accretion correlation. 
This index-mass accretion correlation allows us to evaluate the low limit of black 
hole (BH) mass of compact object in SS 433, Mbht 2 solar masses, using the 
scaling method using BHC GX 339-4 as a reference source. Our estimate of the 
BH mass in SS 433 is consistent with recent BH mass measurement using the 
radial-velocity measurements of the binary system by Hillwig & Gies who find 
that M x = (4.3 ± 0.8) solar masses. This is the smallest BH mass found up to 
now among all BH sources. Moreover, the index saturation effect versus mass 
accretion rate revealed in SS 433, like in a number of other BH candidates, is the 
strong observational evidence for the presence of a BH in SS 433. 
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1. Introduction 

The famous object SS 433 (V 1343 Aql) holds a special place in late twenty century 
Astronomy as the first microquasar discovered in our Galaxy [see reviews by Margon (1984) 
and Fabrika (2004)]. Observations of SS 433 have been carried out in all energy ranges 
during more than 30 years. Its key observational feature is the 162.5 day precession period 
of the jets that is revealed by the lines features. The radial velocity curves of these lines are 
well described by a Kinematical Model which reveals key parameters of the jets v=0.26c, and 
i=79° (Margon 1984). Moreover, Romney et al. (1987) used and combined these results 
with the radio observations of the associated SNR W50 which allowed them to estimate 
distance of 5 kpc to SS 433. 

SS 433 is an X-ray/optical binary with Algol type orbital eclipses. This source is char- 
acterized by two distinct spectral states: the quiescent hard state in which the persistent 
jet flow takes place and the soft state when massive jet blobs are ejected (Fielder et al. 
1987). While the quiescent state has been well studied with numerous X-ray missions, only 
a few massive jet ejection events were seen so far simultaneously with the X-ray soft state 
[Safi-Harb & Kotani (2003), Band et al. (1989)]. Because the ejection of a massive jet blob 
is a rare (in average two times per year) and a short (approximately for ten days) event, we 
only have observations for the part of this jet ejection. 

SS 433 shows many kinds of variability related to regular (orbital and precessional) and 
irregular (flaring) activities. Although precession, binary orbital and nutation periods (162 rf , 
13^.08, 6. d 28 respectively) are well known, but a noticeable variability associated with shorter 
scales is poorly investigated. It is worth noting that the fast variability on time scale of a few 
minutes was investigated by Zwitter et al. (1991) and Goranskij et al. (1987) in the optical 
V-band. More recently X-ray fast variability of 50 seconds time scale during flaring stage 
was found using RXTE observations by Kotani et al. (2002). Power spectrum of SS 433 is 
well presented and approximated by power law (P oc v~ a ) in the range of the 10~ 7 — 10~~ 2 
Hz according to X-ray timing data analysis performed by Revnivtsev et al. (2006). They 
also demonstrated that at frequencies lower than 10 -5 Hz the same variability pattern takes 
place in the optical, radio and X-ray spectral bands (see Fig. 1 of that paper). 

Many questions regarding the complex behavior of SS 433 during outburst states as well 
as the nature of this compact object and its mass are still not answered. However there is 
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no shortage of models which are based on radio, optical and X-ray variations of radiation 
detected from SS 433 [see e.g. Marshall et al. (2002), Fabrika (2004), Safi-Harb & Kotani 
(2003)]. 

The variation in mass estimates of the compact (M x ), optical (M v ) objects and the mass 
ratio (q = M x /M v ) are quite broad. Kawai et al. (1989) and later Antokhina at el. (1992) 
using GINGA observations of SS 433 estimated q ~ 0.15 and q = 0.15 — 0.25 respectively. 
On the other hand Kotani et al. (1996) using ASCA observations found q ~ 0.06 — 0.31 in 
the frame of the precessing jet model with taking into account thermal adiabatic cooling of 
the jets (Brinkmann et al. 1991). Later high-resolution observations by Gics et al. (2002) 
found the presence of absorption lines in the spectrum of the optical A (~A7Ib) supergiant 
companion. These orbital Doppler shifted absorption lines and stationary He II emission 
from the companion allowed to estimate the mass ratio q = 0.35, implying the binary masses 
M x = 4.3±0.8M© and M v = 12.3±3.3M e in SS 433 [see Hillwig & Gies (2008) for details]. 
Thus the average mass ratio inferred from this X-ray data analysis q < 0.25 is smaller than 
that inferred from optical observations q ~ 0.35. 

One comes to the conclusion that in the literature there is a large variation in the mass 
estimates of the compact object ( M x ), secondary star (M v ) and their mass ratio (q = M x /M v ) 
in SS 433. The nature of the compact object was inferred using the mass estimate or its 
upper limit. No other strong arguments were used to determine the nature of the compact 
object in SS 433 which is an eclipsing X-ray binary system, with the primary most likely a 
black hole, or possibly a neutron star [see e.g. Cherepashchuk (2002)]. 

In this work we apply a substantially new approach for diagnosing the nature of compact 
object in SS 433. In §2.1 we present details of radio and X-ray observations of SS 433. 
Analysis of X-ray spectra are shown in §2.2. We discuss X-ray spectral evolution in SS 
433 in §3. X-ray spectral properties as a function of orbital phase are investigated in §4. 
The results of timing and power spectrum analysis are presented in §5. We consider an 
interpretation of observational results and show' our arguments for BH presence in SS 433 in 
§6. We make discussion and concluding remarks in §7. 


2. Observations and Data Reduction 

2.1. Listing of X-ray observations used for data analysis 

We analyzed the archival data collected by PCA / RXTE (Bradt et al. 1993) which 
were obtained in time period from April 1996 to December 2006. These data allow us to 
investigate SS 433 in the broad X-ray energy band (3 150 keV) during quiescent and 



- 4 

outburst states. The RXTE data for SS 433 are available through the HEASARC public 
archive (http://heasarc.gsfc.nasa.gov) at the NASA Goddard Space Flight Center (GSFC). 
As we have already mentioned SS 433 shows continuous (associated with quiet state) and 
sporadic (associated with active state) variability. For investigation of the outburst state and 
for comparing it with the quiescent state we selected only observations during uneclipsed 
intervals of the binary orbital period. In fact, X-ray eclipse occurs around optical primary 
minima at phases |</?| < 0.1. As a result we only used observations taken at interval \ip\ > 0.1 
to exclude the eclipse orbital modulation. In total, this type of observations includes 90 
episodes of outside of eclipse phases. Moreover, 27 observations during eclipses taken at 
different precessional and orbital phases were used for spectral and timing analysis of orbital 
modulation effects. 

Precessional ephemerids were taken from Fabrika et al. (2004). The moment of maximal 
separation between emission lines (T3) was taken to be 73=2443507.47 JD, the precessional 
period -P prec = 162.375 days, the orbital period P or fc=13. 08211 days, the moment of primary 
optical eclipse T o =2450023.62 JD (Goranskij et al. 1998). 

Standard tasks of the HEASOFT/FTOOLS 5.3 software package were utilized for data 
processing. We used methods recommended by RXTE Guest Observer Facility according to 
the RXTE Cookbook (see http://heasarc.gsfc.nasa.gov/docs/xte/recipes/cook book.html). 
For spectral analysis we used PCA Standard 2 mode data, collected in the 3 - 20 keV energy 
range. The standard dead time correction procedure has been applied to the data. To 
construct broad-band spectra, HEXTE data have been also used. We subtracted background 
corrected in off-source observations. To exclude the channels with largest uncertainties only 
data in 20 - 150 keV energy range were used for the spectral analysis. In Table 1 we list 
the groups of RXTE observations covering the source evolution from quiescent to outburst 
states. We also used public data from the All-Sky Monitor (ASM) on-board RXTE. The 
ASM light curves (2-12 keV energy range ) were retrieved from the public RXTE / ASM 
archive at HEASARC l . 

In the present Paper, we have analyzed X-ray spectra during quiescent and outburst 
states with reference to simultaneous radio and optical observations. The monitoring RAT AN- 
600 Radio Telescope (2-8 GHz) data in the 1996 - 2006 period w r ere available through the 
public archive (http://cats.sao.ru/~satr/BH). We also used radio observations by Green 
Bank Interferometer, NRAO 2 obtained from 1996 to 1998 at 2.25 and 8.3 GHz and si- 
multaneous V-band photoelectric photometric observations. Details of optical telescopes, 


1 http: / /xte. mit.edu/ASMJc.html 

“http: / / www . gh . nrao . edu/ fdocs / gbi/arcgbi 
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reduction techniques and compilation methods are given by Goranskij et al. (1998). 

Additionally, we analyzed the INTEGRAL/IBIS/ISGRI spectra in flaring state (2004) 
of SS 433, that were coordinated with the RXTE observations. We have used the version 
8.0 of the Offline Science Analysis (OSA) software distributed by INTEGRAL Science Data 
Center [ISDC, http://isdc.unige.ch, Corvoisier et al. (2003)]. 

We also present a comparison of the SS 433 data with that for GRS 1915+105 obtained 
during BeppoSAX observations. We used two BeppoSAX detectors (Medium- Energy Con- 
centrator Spectrometer (MEGS) and a Phoswich Detection System (PDS)) for this analysis. 
The SAXDAS data package was utilized for performing data analysis. We process the spec- 
tral analysis in the good response energy range taking into account satisfactory statistics of 
the source: 1.8 10 kev for MECS and 15 150 keV for PDS. 


2.2. Spectral analysis 

SS 433 has long been of great interests in X-ray Astrophysics, and has been observed 
early on with many satellites such as HEAO-1 (Marshall et al. 1979), EXOSAT (Watson et 
al. 1986), Tenma (Matsuoka et al. 1986), and Ginga (Kawai et al. 1989). Using HEAO-1, 
Marshall et al. were the first to demonstrate that SS 433 is an X-ray source. The HEAO-1 
continuum was sufficiently modeled as thermal bremsstrahlung with kT = 14.3 keV, and 
emission due to Fe-K was detected near 7 keV. The ASCA satellite, which carried X-ray 
CCD cameras for the first time, detected many pairs of Doppler-shifted emission lines from 
ionized metals, such as Si, S, Ar, Ca, Fe, and Ni, originating from the twin jets (Kotani et al. 
1994). The emission lines were also resolved with the Chandra HETGS, which were found to 
have Doppler widths of 1000-5000 km s' 1 (Marshall et al. 2002; Namiki et al. 2003; Lopez 
et al. 2006). The broad band continuum (up to 100 keV) is approximated by a thermal 
bremsstrahlung spectrum with a temperature of 10-30 keV, depending on whether SS 433 
is in or out of eclipse (Kawai et al. 1989; Cherepashchuk et al. 2005). Additional complex 
features were detected from the XMM spectra, however, which could be Compton-scattered 
emission from the jet base (Brinkmann et al. 2005) or an iron-K absorption edge due to 
partial covering (Kubota et al. 2007). From the width of an eclipse in the 25-50 keV band 
with INTEGRAL, Cherepashchuk et al. (2007) and Krivosheyev et al. (2009) propose that a 
hot extended corona around the accretion disk is responsible for the hard X-ray emission via 
thermal Comptonization with a temperature of 20 keV. High-quality X-ray spectra covering 
the broad band are critical in establishing an interpretation of the high energy spectra of SS 
433. 



In our study we model the broad band source spectra in XSPEC using an additive model 
consisting of sum of the so called bulk motion Comptonization ( BMC) model and two Gaus- 
sian line components. The BMC model is a generic Comptonization model which can be 
applied to upscattering of soft photons injected in a hot cloud. This model consists of two 
parts where the first part is a direct blackbody (BB) component and the second one is a 
convolution of the fraction of the BB component with a broken power law, the upscattering 
Green function. The spectral index of the blue wing a is much smaller than that of the red 
wing a + 3. The shape of the Green function (broken power law) is generic and independent 
of the type of Comptonization, thermal or nonthermal. The name of the model (BMC) has 
only a historical sense dating to 1997 [see Titarchuk et al. (1997)] when the model was hrst 
applied to case of the bulk motion Comptonization. However this model can be applied to 
any type of Comptonization, thermal or nonthermal but it should be, in principle, combined 
with exponential cutoff which is related to average plasma energy, for example, plasma tem- 
perature for the thermal Comptonization kT e or kinetic energy of the matter in the case of 
the converging (bulk inflow) Comptonization. In this Paper we consider a scenario related 
to our model (see Fig. 1) where the Compton cloud along with converging flow are located 
in the innermost part of the source and a Keplerian disk extends from the Compton cloud 
(CC) to the optical companion. 

As we point out ASCA and Chandra detected many lines of various elements in the soft 
X-ray band of spectrum of SS 433. Particularly, iron lines, Fe XXV - Fe XXVI dominate 
at energies 6.5<E<7 keV and these lines show a double structure due to jet Doppler shifts. 
In addition to iron line emission one can see the line emission related to hydrogen - and 
helium-like ions of Mg, Si, S, Ar, Ca and Ni which display a double structure also. These 
line signatures indicate that the lines are formed in the relativistic jet configuration. Along 
with these lines there is an appreciable emission feature at 6.4 keV which is visible in the 
X-ray spectrum of SS 433 [Kotani et al. (1996); Seifina (2000)]. This line is not subjected 
to Doppler shifting. Thus we want to emphasize that using the forms of these lines we can 
see features of “moving” and “stationary” structures of material surrounding SS 433. 

However the identification and precise theoretical reproduction of the line composition 
with RXTE is problematic because of its low energy resolution. As a test trial we added 
one Gaussian component to fit the spectrum varying the width and normalization of the line 
and we found that the width a of this Gaussian feature roughly ranges from 0.3 to 1 keV. In 
quite a few cases the spectral fits using one Gaussian component provide very wide residuals 
extended from 6 to 9 keV. However, after adding a second narrow Gaussian component (in the 
6-9 keV range) the fit quality has been significantly improved. The energies of the first and 
second Gaussian components Ey ne i and Ey ne 2 are presented in Table 3. For the first Gaussian 
Eu ne 1 changes from 6.5 to 6.9 keV while the range of the second Gaussian varies from 7.1 
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to 9 keV. In some cases we see a wide residual taking place around 20 keV which can be a 
signature of ’’high temperature bbody-like” spectral component of temperature in the range 
of 4-5 keV. Thus we use our XSPEC model as wabs*(bmc+Gaussian+Gaussian+bbody ) for 
fitting of SS 433 spectra. In particularly we use a value of hydrogen column N # = 1.2 x 10 23 
cm 2 which was found by Filippova et al. (2006) in calculations of XSPEC model wabs. 

The best-fit parameters of the source spectrum are presented in Tables 3 6. For the 

BMC model the parameters are spectral index a (photon index F = a + 1), color tempera- 
ture of the blackbody-like injected photons kT, log(A) related to the Comptonized fraction 
/ [/ = A/ (l + A)] and normalization of the blackbody-like component Nf mc . We find that 
color temperature kT is about 1 keV for all available RXTE data and thus we fix a value of 
kT at 1 keV. When the parameter log(A) 3> 1 we fix log(A) = 2 (see Table 3, 5), because the 
Comptonized fraction f = A/( 1 + A) — >1. The variations of A do not improve fit quality any 
more. A systematic error of 1% has been applied to the analyzed X-ray spectra, we applied 
this systematic error to the analyzed RXTE spectra in accordance to the current version of 
RXTE Cookbook (see http://heasarc.gsfc.nasa.gov/docs/xte/recipes/cook book.html) fol- 
lowing to recommended methods by RXTE Guest Observer Facility. 

In Figure 1 we present a suggested geometry of X-ray source in SS 433 (see more 
explanation of geometry below). 


Similar to the ordinary bbody XSPEC model, the normalization Nb mc is a ratio of the 
source (disk) luminosity to the square of the distance D 


N, 


bmc 


(lO 39 erg/s) ( 


lOkpc 

D 


( 1 ) 


The adopted spectral model shows a very good performance throughout all of the data 
set used in our analysis. Namely, the value of reduced y 2 -statistic Xred = X 2 /^dof > where 
Ndof is a number of degree of freedom for a fit, is less or around 1.0 for most observations. 
For a small fraction (less than 2%) the spectra with high counting statistic xf. cd reaches 1.5. 
However, it never exceeds a rejection limit of 2.0. 


3. Evolution of SS 433 X-ray spectra during outbursts 

We have established general tendencies of quiescent-outburst behavior of SS 433 based on 
spectral parameter evolution of X-ray emission in energy range (3-150) keV using RXTE/P- 
CA/HEXTE data. We have also found correlation between X-ray emission of SS 433 and 
radio patterns observed by RATAN-600 in the range 1-11 GHz and by Green Bank Interfer- 
ometer at 2.25, 8.3 GHz. 



8 - 

We identified common features of ouburst behavior of SS 433 based on three available 
outbursts and on two outburst decay sets. According to general BH state classification [see 
for example Shaposhnikov k Titarchuk (2009), hereafter ST09] SS 433 is mostly seen in the 
intermediate (IS) state. Different patterns of SS 433 X-ray spectra are listed at Table 3: i) 
IS spectra with T ~ 1.9 — 2.2 (see also lower panel of Fig. 2); ii) sum of BMC spectrum with 
T ~ 2.2 — 2.3 and high temperature “bbody” component found in 10-50 keV energy range 
(see Figs. 3-4). 

Ten days before radio outburst, SS 433 being in IS a stable low soft (ASM) X-ray flux 
(see Fig. 5) is followed by a X-ray flux rise reaching its maximum just two days before the 
radio flare. Probably at the moment of radio peak, MJD=50890 the X-ray flux reaches its 
minimum (BMC normalization, Nbmc ~ 0.7 x 10~ 3 ) and when X-ray spectrum becomes 
harder (photon index T ~ 1.9). It is worth noting that spectral index of radio emission 
otradio has its maximal value at a time close to X-ray outburst (see Fig. 5, at MJD=50886). 
Note the radio spectrum is harder just a few days before the radio maximum. Then, after 
several days of the radio maximum, we see the X-ray flux increases while the radio flux 
substantially decreases, at least by factor 2 (see Fig. 5 and Table 3). 

In Figs. 6, 7 we demonstrate an evolution of the radio, X-ray flux and X-ray spectral 
parameters Nb mc and F with time for different time intervals. Notably in Fig. 6 we present 
the best-fit color temperature related to “high temperature BB” bump which is sometimes 
seen in IS spectra. Also in Figs. 5—7 we show a transition between intermediate and soft 
states given that photon index varies between 1.9 and 2.4. 

Thus in Figures 5—7 we present the behavior of SS433 during corresponding outbursts 
and spectral transitions using i. the BMC normalization, ii. the photon index, iii. ASM 
X-ray flux, and sometimes when they are available, iv. GBI/RATAN-600 radio flux, v. radio 
spectral index (in case of GBI observations) and the color temperature (in the case of the 
occurrence of ‘’high temperature Black body” -like feature). 

The diagram of photon index T versus BMC-normalization forms the track shown in 
Fig. 8 for all observations except ones during orbital eclipses. The rising part of the correla- 
tion is seen in both precession periods and outburst transitions. The upper saturation part 
of the correlations is only seen during outburst transitions. Note that in the cases of different 
outbursts of SS 433 the levels of photon index saturation are also different. For example, 
the 2004 outburst (blue square) shows the saturation level around F ~ 2.3, while the 2005 - 
2006 outbursts (crimson and red circles) saturates at lower level of F ~ 2.2. Comptonization 
fraction / shown in the right hand panel of Fig. 8 is high in most of cases. This means that 
in most cases the soft disk radiation of SS 433 is subjected to reprocessing in a Compton 
Cloud and only a small fraction of disk emission component (1 — /) is directly seen. Thus 
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the energy spectrum of SS 433 during all states is dominated by a Comptonized component 
seen as a power-law hard emission in the energy range from 3 to 70 keV, while the direct 
disk emission is not seen in all detected spectra. 


4. Spectral properties as a function of orbital phase 
4.1. Detection of “high temperature BB-like” bump 

The “BB-like” feature has been found in 24 spectra of SS 433 among of all available data 
(R2, R3, R5-R.7). Note that during 2004 outburst (R5 set) “high temperature BB-like” bump 
found in 10-40 keV energy range of RXTE has been also detected by ISGRI/IBIS detector 
onboard INTEGRAL satellite during simultaneous with RXTE observation (ID 90401-01- 
01-01). For two time intervals (MJD=52225-52238, 53579-53588) this feature is detected in 
spectra near the primary eclipse. It has been visible before and after the primary eclipse, 
while it is not seen during the central phases of the primary eclipse. This screening effect 
with respect to orbital phases is clearly seen in two lower panels (E3 and E4) of Fig. 9. Blue 
vertical strip marks the interval of the primary eclipse made according to optical ephemerids 
(Goranskij et al. 1998). Points marked with rose oreol correspond to spectra fitted by 
the model which includes ’’high- temperature BB-like” component (for details see Fig. 3, 4 
and Tables 4-7). Red points in Fig. 9 ( see two lower panels ) correspond to observations 
during primary eclipse when the ” BB-like” component is not seen, while it has been found 
before and after the eclipse. Although observations are distributed more or less uniformly 
with precession phase, the ”BB” bump feature is seen generally at ’’disk face-on” position 
[ip = 0.75 (E3), ip = 0 (E4), Fig. 9], when the precessing disk is open and the innermost part 
of the disk is more visible by the Earth observer. 

Thus we argue that the ” BB-like” feature is better detected at ’’disk face-on” position 
(0 — *• 0) than that during other precession phases. Moreover we have found the total eclipse 
of ”BB” bump during the primary eclipse when the compact object is obscured by the normal 
star. In other words, we suggest that localization of ”BB” bump feature emission should 
be better identified during the direct observation of the innermost source region near the 
compact object (CC). 

However this “BB” component can be smeared out by photon scattering if optical depth 
of Compton cloud is much greater than 1. That could be the case during observations for 
MJD=50191-50194 (El) and 50897-50907 (E2) (see Fig. 9). 



4.2. Softening of spectra at some particular orbital phases 


The relative softening of the spectrum was found in the range of ip = 0.95 — 1.1 for two 
orbital cycles of SS 433 (see Fig. 9). Softening of X-ray spectra is seen when the compact 
hard emitting region, presumably related to the central parts of the disc, obscured by the 
donor (optical) star and thus the normalization of X-ray spectrum is suppressed due to the 
eclipse. This relative spectral softening is clearly seen at M JD=50 191-50194 (El), 50897- 
50907 (E2) (see two upper panels of Fig. 9) when the noticeable increase of photon index at 
MJD=50192.7 (T = 2.43 ±0.03) occurs with respect of that during the rest of eclipse interval 
(r = 2.11 ± 0.01) (see El panel in Fig. 9). The similar behavior of index is seen in E2 panel 
where the noticable increase of photon index T = 2.39 ± 0.01 takes place at MJD=50899.7 
with respect of that when F = 2.13 ± 0.01 during the rest of the eclipse interval. 


5. Timing analysis 

The RXTE light curves were analyzed using the powspec task from FTOOLS 5.1. The 
timing analysis RXTE / PCA data was performed in the 2-15 keV energy range using the 
binned mode. We generated power density spectra (PDS) with 16-second time resolution. We 
subtracted the contribution due to Poissonian statistics. To model PDS we used QPD/PLT 
plotting package. 

We also analyzed optical and radio light curves to compare them with X-ray curve. 
Timing analysis of all available data of SS 433 demonstrates power-law PDSs of index around 
1.5 in three energy bands (see Fig. 10): X-ray (blue), 2.25 GHz radio (black), V-optieal (red). 
This result is in an agreement with previous results by Revnivtsev et al. (2006) who found 
that the SS 433 power spectra of radio, optical and X-ray variabilities continue with the 
same power law from 10~ 7 Hz up to 10 '' Hz. 

Such art extended behaviour of power spectra is well known for a number of X-ray 
binaries. For instance, Churazov et al. (2001) found that Cyg X-l demonstrated the similar 
power law component at 10~ 5 — 10 Hz in X-ray energy band. An extensive study of low mass 
X-ray binary systems showed that such a power law behavior is common and often observed 
in LMXBs [see e.g. Gilfanov & Arefiev (2005) and Titarchuk et al. (2007)]. 

The observed variability of SS 433 in different spectral bands was interpreted in the 
framework of self-similar accretion disk variations suggested by Lyubarskii (1997). In par- 
ticular, he proposed that when variable mass accretion rate generates energy release in the 
disk far away from a compact object then the corresponding radiation should lie seen in 
optical and UV energy range. Whereas the energy emitted in X-ray should be released near 
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the compact object. In SS 433 presumably these mass accretion rate fluctuations are de- 
tected and observed in power density spectrum. In which one can also see that the X-ray 
variability correlates with the optical and radio emissions (see Fig. 10). 

A study of orbital modulation allow 11 s to detect some changes of power spectrum of 
SS 433 in the energy range 3 15 keV during the primary eclipse. In the lower panel of 

Fig. 11 we show details of the typical evolution of X-ray timing and spectral characteristics 
during the primary eplipse. The top panels of this Figure demonstrate V-band optical light 
curve, BMC normalization and photon index T for R6 set as a function of time (see Table 1 
for details of these observations). In the plot of photon index T versus time points A and B 
correspond to moments MJD=53581 and 53585, during and after eclipse respectively. Blue 
strips mark eclipsed intervals of light curve. It is worth noting that the optical and X-ray 
fluxes significantly drop during the primary eclipse (see two upper panels and Table 3). 

In Fig. 11 PDSs ( left bottom panel ) are plotted along with energy spectrum E x F(E) 
( right bottom panel) for two observational points A ( 91103-01-03-00 , blue) and B (91103- 
01-07-00, red). The differences of the presented PDS and energy diagrams related to events 
A and B are noticeable. PDS during eclipse related to point A (blue histogram) indicates 
some weakening of power at 0.001 0.01 Hz in comparison with that for PDS outside of 

eclipse, at point B (red histogram). The corresponding spectrum during eclipse (blue line) 
demonstrates a significant decrease of total flux and also of the BMC normalization, at least 
by factor 2. 


6. Interpretation of observational results 

Before to proceed with the interpretation of the observations let us to briefly summarize 
our results as follows: i. The spectral data of SS 433 are well fit by a BMC model plus 
two Gaussian and “high temperature bbody” components (see §2.2). ii. Power law spectral 
index of BMC component rises and saturates with an increase of the BMC normalization 
(disk flux). The photon index saturation levels of the 2004 and 2005 -- 2006 outbursts are 
about 2.2 and 2.3 respectively (see Fig. 8). iii. There is a total eclipse of ’’high temperature 
BB-like” feature in the photon spectrum because this spectral feature is presumably formed 
in the innermost part of the X-ray source and thus it should not be seen during the primary 
eclipse, iv. The spectral index (hardness) clearly changes during the primary eclipse (see 
panels El, E2 in Fig. 9 and one of upper panels in Fig. 11). This can be explained if there 
is a slight variation in spectral shape with radius through X-ray spectral formation region. 



6.1. Index saturation as a BH signature. BH mass in SS 433 


We have firmly established that the photon index of the resulting spectrum of SS 433 
saturates with the BMC normalization N hmc oc L d or with disk mass accretion rate. ST09 
presented strong arguments that this index saturation is a signature of converging flow into 
BH. In fact, the spectral index a is inverse proportional to Comptonization parameter Y 
which is proportional to product of average number of photon scattering in Compton cloud 
N sc and mean efficiency of gaining energy at any scattering rj. Because the index saturates, 
it means rj is an inverse of N sc . But in general for an optically thick Compton cloud of 
optical depth r » 1 average number of scatterings N sc is proportional to r 2 or r [see e.g. 
Laurent & Titarchuk (1999)]. In the case of thermal Comptonization A T iSC oc r 2 and then r] 
should be inverse proportional to r 2 . But in the thermal Comptonization case the efficiency 
r] is proportional to electron temperature kT e only [see e.g. Sunyaev & Titarchuk (1980)] 
and it is not clear why there could be such a sophisticated tuning in the Compton cloud 
that kT e steadily decreases with r as r~ 2 . On the other hand N sc oc r and 77 oc 1/r in 
the converging flow case [see e.g. Titarchuk et al. (1997) and ST09] and thus the index 
saturation is a natural consequence of these laws for N sc and 77 in the case of converging flow 
[see Titarchuk & Zannias (1998) for the solution of the full relativistic kinetic equation]. 
Furthermore ST09 demonstrate that for a number of BHC sources index saturates not only 
when the mass accretion rate increases but also when QPO frequency increases. These 
index correlation and saturation with QPO frequency are a direct evidence that the index 
saturation takes place when the size of Compton cloud (CC) becomes more compact (namely 
the index saturates when compactness of the emission area increases). 

Titarchuk et al. (1998) demonstrated that the CC compactness increases with mass 
accretion rate, M. Moreover in the BH case the shape of the emergent spectrum saturates, 
when m increases, to the spectra for which the power-law indices F stay almost the same 
and these F are always below 3 (see ST09 and Fig. 8 in this Paper). In fact, the asymptotic 
values of F in SS 433 vary in the ranges ~ 2.2 - 2.3 (see Fig. 8). Thus, we argue that 
the SS 433 X-ray observations reveal the index saturation vs mass accretion which can be 
a sign of a converging flow in this source indicating a BH presence. This index vs BMC 
normalization (proportional to disk mass accretion rate) correlation can, in principle, allows 
us to estimate a lower limit of BH mass in SS 433. 

To do BH mass estimate in SS 433 we choose to scale index-normalization previously 
found in GX 339-4 (see ST09) with that in SS 433 (see Fig. 12). We can proceed with 
this scaling if these two correlations are self-similar. In other words when given correlations 
have the same the index saturation levels and slopes as a function of normalization which is 
proportional to mass accretion rate in the disk. 
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Then the value of the scaling coefficient provides us BH mass estimate. Note that the 
index-normalization correlation curve for 2003 decay data of GX 339-4 (taken from ST09) 
is self-similar with that we find for SS 433. 


Use of scaling method for these two correlations provides the scaling coefficient S n = 
14.9 ± 0.7 which is a ratio (see ST09): 


s n — 


N r _ M r d‘f 
Ni. ~ JI^P r fG 


( 2 ) 


where M r is BH mass of the reference source (GX 339-4), M t is BH mass of the target 
source (SS 433), geometric factor fc = (cos9) r /(cos0) t , inclination angles 9 r , 9 t and d r , d t 
are distances to the reference and target sources. Here the subscripts r and t are related to 
the reference and target sources respectively. One can, in principle, estimate BH mass M t in 
SS 433 using the relation (2). However a value of geometric factor fc = (cos8) r /(cos9) t is 
not obvious because cos 9 r which is an inclination angle of X-ray emission area in GX 339-4 
is not known. Also it is not clear if we can use the inclination of the whole disk in order 
to estimate the inclination of the innermost part of the disk which size is about 4 orders of 
magnitude less than that of the disk. This is one of the main problems to correctly evaluate 
.fa- 


in order to estimate (cos#) r = (cos 9)gx we can apply the same formula (2) but using 
the target and reference sources GX 339-4 and GRO J1655-40 respectively. We take values 
of M r , M t , d r , d t and cos 9qro from ST09 (Table 5 there) and then we obtain the lowest 


limit of 


(cos 0) G J>C> 0.4 


(m r /6.3)((ij/5.75 kpc) 2 [(cosfl) G Ko/ cos70°] 
(m t /12.3)(d r /3.2 kpc) 2 


( 3 ) 


Now we can resolve formula (2) with respect to m t = m ss where we use GX 399-4 and SS 
433 as the reference and target sources respectively. Thus we obtain 


> (A/5.0 kpc) 2 [(cosfl) o , Y /0.4](m r /12.3) 

(D r j 5.75 kpc) 2 [(cos#)., s /0.17] (sw/14.9) 


( 4 ) 


where we use a value of cos 9 SS ss 0.17 estimated by Cherepashchuk et al. (2005). 

This is the lower limit mass estimate for SS 433 found using the ST09 scaling method. 
This estimate is consistent with the BH mass value Mss = (4.3 ± O.8)M 0 obtained using 
optical observations of SS 433 [see Hillwig &; Gies (2008), hereafter HG08]. BH mass 
estimates for the reference (GX 339-4) and target (SS 433) sources are summarized in Table 7. 



6.2. Upper limit of the size of “High temperature BB-like” bump region 


Now wc try to estimate the upper limit of the size of the region A Ibb from where the ” 
high temperature BB-like” photons originates. We consider the screening of the ” BB-like” 
bump emitting region by the optical star with radius R st . 


Based on BII mass estimate m x by HG08, one can determine the mass of donor (optical) 
star using the mass function of normal star f v (m ) for SS 433. Namely, 


fv(m) 


m x 3 sin 3 i 
( m x + rn v ) 2 


= 1.038 x 10” 7 A' w 3 * P orfc (l - e 2 ) 3/2 , 


(5) 


where K v = 58.2 ±3.1 km s" 1 is semi-amphlitude of radial velocity, P OT ), = lS'bOS is orbital 
period in days, e is eccentricity (e ~ 0), i = 78°. 8 is inclination angle of system [see Margon 
(1984)]. Thus mass of the donor star rn v can be obtained as 


m, 


I m% sin 3 i 
fv(m) 


rn x ~ 12.3 M® 


( 6 ) 


if one resolves Eq.(5) with respect to m v and uses the BH mass values found by HG08 (see 
above). Then the resulting mass ratio is q = m x /rn v < 0.35 [see HG08 and Filippova et al. 
(2006)]. 

Using the donor star mass one can calculate the binary separation and Roche lobe 
radius. The resulting binary separation from Kepler’s third law is a ~ 6OP, 0 . Because the 
donor (optical) star in SS 433 fills its Roche lobe R st = R cr , the Roche lobe radius of the 
donor star can be found (compare with Eggleton, 1983) 


Rgt ~ a x 


0 . 49 i /~ 2 / 3 

0.6g -2 / 3 + ln(l + q~ U 3 ) 




28 R e 


(7) 


But the upper limit of the ” BB-like” bump forming region can be estimated as 

A Ibb ~ - — p — — — 2 * R, st < 6.8 Rq. (8) 

^orb 

using that ” BB-like” bump forming region is screened, during A t BB < 2.5 days, by a donor 
(optical) star of radius R st and rotating with the orbital velocity v = 2ira/ P or i, (see Figs. 9, 
13). It is worth noting that A Ibb, the real size of u BB-like” component emission, should be 
significantly smaller than that we estimate by inequality (8). Here we have just found the 
upper limit of A Ibb based on the direct observational estimates. 
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7. Discussion and Summary 

We derive a lower limit of 2 solar masses for the compact object in SS 4333 which clearly 
does not rule out a neutron star - this limit alone does not give strong observational evidence 
for a black hole. This is the lower limit of the central object mass is obtained using sliding 
the index-mass accretion correlation of SS 433 with respect to that for GX 339-4. In fact, this 
is an example that the sliding coefficient provides the lower limit of the mass of the central 
object. We remind a reader that this correlation of photon index and mass accretion rate 
was obtained using the so called BMC model of XSPEC which is the first principle radiative 
transfer model. Namely we consider and confirm applying the XSPEC fitting procedure that 
the emergent spectrum of the Compton cloud (CC, the innermost part of a BH) is a result 
of upscattering of soft, presumably disk, photons off CC energetic electrons. 

The crucial characteristic of the correlation is the saturation plateau which is a signa- 
ture of the converging flow into central object because the spectral index is an inverse of 
Comptonization parameter Y which is a product of number of scatterings N sc and average 
fractional energy change per scattering 77 (see also details of explanation of index saturation 
in §6.1). Only for converging flow N sc is proportional to optical depth r of the cloud and 
rj is inverse proportional to r when r is greater than one. Thus for large r (mass accretion 
rate) it is natural to expect the index saturation as a signature of a BH and in fact, we really 
see this saturation in the case of SS 433. In the other words the main evidence for black 
hole is the index saturation with mass accretion rate. In contrast, in neutron star sources 
the photon index is almost constant, concentrating near 2, i.e, it does not change when mass 
accretion rate increases [see Farinelli k Titarchuk (2010)]. 

Moreover Hillwig k Gies (2008) found, using optical observations of SS 433, that the 
mass of central object is 4.3 ± 0.8 solar masses which can be the case only for a black hole. 
It is also worth noting a recent paper by Kubota et al. (2010) who argue these mass values 
should be considered as upper limits and taking into account the heating of the donor star 
may lower the derived masses to 1.9M® < M x < 4.9 M§. 

We find that the BMC model along with an additional “high temperature bbody-like” 
component allow us to consistently describe X-ray spectra of SS 433 and derive physical con- 
ditions during all states of SS 433. Moreover, our approach for data fitting demonstrates that 
the SS 433 spectra becomes softer and finally photon index saturates when mass accretion 
rate increases during X-ray outburst. This index behavior is in agreement with that ST09 
and TS09 previously find in a number of other BHC sources. Furthermore, the fits, using our 
spectral model, allow us to demonstrate that in some cases an additional “high temperature 
bbody -like" component of color temperature of 4-5 keV appears in the data (see Figs. 3-7 
and Tables 4-6). This “bbody-like” feature arises only during outburst transition of SS 433 
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and it presumably originates in the innermost part of accretion flow (see discussion in section 

6 . 2 ). 

Recently we found a signature of this "high temperature bbody-like” bump in RXTE 
spectral data for BHC GRS 1915+105 (see TS09). Also we find that the same kind of 
“blackbody-like” bump is seen in BeppoSAX observations of the intermediate state of GRS 
1915+105 (April 21, 2000, ID=209850011). The data along with the model are presented 
in Figure 14 (top panel). Consequently the presence of the additional ‘‘BB-like” bump 
spectral feature is not an instrumental effect of RXTE but it can be a real observational 
feature of the observed spectra of BHCs. Furthermore, we also detected the same ‘bbody- 
like high temperature” feature in the spectrum of SS 433 using IBIS/ISGRI detector onboard 
INTEGRAL satellite during simultaneous observations with RXTE (ID=90401-01-01-01, R.5 
set) (see Fig. 14 bottom panel). Thus it is quite possible that this “high temperature bbody- 
like ” feature can be the universal signature of BH spectra. 

Also one can argue that the presence of “high temperature bbody-like” feature may be 
an artifact of having chosen the BMC model. In fact, the presence or absence of the bump 
can be related to the model. But the BMC model is a first principle Comptonization model 
and therefore any residual using this model can be considered as an observational signature 
of other component presented in the spectrum. Moreover Laurent & Titarchuk (1999) and 
Laurent & Titarchuk (2010) demonstrate using Monte Carlo simulation that X-ray emergent 
spectra of a black hole for all spectral states are well fitted by BMC model if the effect of pair 
production is not taken into account. But in the spectra of the intermediate and soft states 
this additional “BB features appears if the non-linear effect of pair creations in very vicinity 
of a black hole horizon is taken into account [Laurent & Titarchuk (2010) in preparation]. 

One can say that such a bump seems very similar to the “Compton-bump” expected by 
reflection models in black hole systems (both AGN and stellar-mass objects), e.g. see Miller 
(2007). In fact, in our previous paper TS09 (see section 4.2 there) we discuss this issue in 
detail. Here we should once again point out that in the observations this “high temperature 
BB-like bump appears only when the photon index is higher than 2. On the other hand 
Laurent & Titarchuk (2007) demonstrate using Monte Carlo simulation and theoretical 
arguments that the reflection bump never appears in the emergent spectra if photon index 
of the spectrum F is higher than 2. Namely, there are not enough photons in the incident 
spectrum at high energies, if F > 2, to be reemitted into lower energies due to scattering and 
recoil effects. In fact, as one can see from Tables 4-6, that in all spectra where we detect this 
“high temperature BB-like” (~ 20 keV) feature index of the hard BMC component F > 2 
and thus the appearance of this bump in the observed spectra cannot be explained by so 
called “reflection” effect. 
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TS09 argue that this “high temperature BB-like” bump can be a result of the grav- 
itationally redshifted annihilation line which is originally formed in the very vicinity of a 
BH horizon. This feature (bump) should be seen in the IS and HSS spectra only because 
in these states mass accretion rate is high enough to provide conditions to form the strong 
annihilation line and also to observe this feature through accreting material of relatively low 
plasma temperature (see more details in TS09, §4.2 there). 

It is worth noting that a time delay of the radio peak with respect to that in X-rays 
is two days for outburst rise (Revnivtsev et al. 2006). We find that during this X-ray 
flare, and then in radio one, the object transits from IS to LHS, which is quite unusual for 
spectral transitions of BHC. This particular behavior of SS 433 is different from other BH 
sources probably because the jet is the dominant emission component. In general, when a 
BH goes to outburst it leaves a quiescent state, enters to LHS and then it goes to IS-HSS. 
However, the microquasar SS 433 is the only BHC which stays, most of the time in IS with 
a rare short transition to LHS accompanied by the radio flare. A delay of radio peak with 
respect to X-ray peak is known during outburst rise transition in many X-ray BHC. However 
a time delay value is different for each of the sources. For instance, it is about 20 days for 
GRS 1915+105 and about 2 days during outburst in SS 433. If in the GRS 1915+105 case 
the index saturation value is about 3, whereas in the case of SS 433 the index saturates to 
F ~ 2.3. Moreover the direct soft component, which is usually associated with the disk, is 
clearly seen in the GRS 1915 soft spectra and it is weak or absent at all (at the level of 
detection) in most of SS 433 spectra (see e.g. Table 3). 

As conclusions we formulate the following. We analyze the state transition data from 
SS 433 collected using RXTE observations. We examine the correlation between the photon 
index of the Comptonized spectral component and its normalization which is presumably 
proportional to disk mass accretion rate (see Fig. 8). We find that broad-band energy spectra 
during all spectral states are well fit by the XSPEC BMC model for the continuum and by two 
(broad and narrow) Gaussian line components. In addition to these model components we 
also find a strong feature of ”blackbody-like” bump which color temperature is in the range 
of 4-5 keV in 24 intermediate state (IS) spectra of SS 433. 

Furthermore the application of our spectral model to the SS 433 data allows us to 
establish the saturation of the photon index vs BMC normalization, which scales with disk 
mass accretion rate at value around 2.3 (see Fig. 8). In addition, an application of the scaling 
method (see ST09) allows us to estimate a lower limit of compact object mass in SS 433 
{M x > 1.8M 0 ). 

A high value of the Comptonized emission fraction / obtained, in the framework of 
BMC model, gives us a strong evidence of significant reprocessing of X-ray emission of the 
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disk photons in SS 433 which is also in agreement with power density spectra which reveal 
a pure power law (so called “red noise”). 

It is important to emphasize that the index saturation effect in SS 433 now seen, is a 
BH signature and was recently found in a number of BHC sources (ST09). Moreover the 
detection of the so called “high temperature BB-like” bump” (which could be a gravitation- 
ally redshifted annihilation line) in SS 433 is also found in GRS 1915+105 data by different 
space missions ( RXTE , BeppoSAX and INTEGRAL). 

We acknowledge Vitalij Goranskij and Tatyana Irsmambetova who kindly provide us 
optical data. We also thank RXTE , INTEGR AL team and also Sergej Trushkin for providing 
us X-ray and radio data correspondingly. We appreciate comments and thorough editing of 
the paper made by Chris Shrader. We are very grateful the referee for his/her valuable 
comments and corrections of the content of the paper. 


REFERENCES 

Antokhina, E. A., Seifina, E. V., & Cherepashchuk, A. M. 1992, Soviet Astronomy, 36, 143 
Band W. et al., 1989, ApJ, 336, 937 

Bradt, H.V., Rothschild, R.E. & Swank, J.H. 1993, A&AS, 97, 355 

Brinkmann, W., Kawai, N., Matsuoka, M., Fink, H. H., 1991, A&A, 241, 112 

Cherepashchuk, A.M. et al. 2005, A&A, 437, 561 

Cherepashchuk, A.M. 2002, Space Science Reviews, 102, 1 

Churazov, E., Gilfanov, M. & Rcvnivtsev, M. 2001, MNRAS, 321, 759 

Corvoisier, T.J.-L. et al. 2003, A&A, 411, L53 

Eggleton, P.P., 1983, ApJ, 268, 368 

Gies, D. R., McSwain, M. V., Riddle, R. L., Wang, Z., Witta, P. J. & Wingert D. W. 2002, 
ApJ, 566, 1069 

Gilfanov, M. & Arefiev, V., 2005, astro-ph/0501215 

Goranskij, V.P., Kopylov, I.M., Rakhimov, V.Yu. et al., 1987, Cornmun. Spec. Astrophys. 
Obs. 52. 5 



19 

Goranskij V.P. et al., 1998, Astron. Rep., 42, 209 
Fabrika, S.N. 2004, Astrophys. Space Pliys.Rev., 12, 1 
Farinelli, R. k Titarchuk, L. 2010, submitted to AkA 

Filippova, E., Revuivtsev, M. Fabrika, S., Postnov, K. k Seifina, E. 2006, AkA, 460, 125 

Fielder, S. et al. 1987, ApJ, 94, 1244 

Hillwig, T. C. k Gies, D. R. 2008, ApJ, 676, L37 

Hynes, R.. I., Steeghs, D., Casares, J., Charles, P. A., & O’Brien, K. 2004, ApJ, 609, 317 

Kawai, N., Matsuoka, M., Pan, H., k Stewart, G. C. 1989, PASJ, 41, 491 

Kotani, T., Band, D., Denissyuk, E.K., Kawai, N., Kinugasa, K., Namiki, M., Safi-Harb, S., 
k Trushkin, S. 2002, in ASP Conf. Ser. 279, Exotic Stars as Challenges to Evolution, 
ed. CA Tout k W. Van Hamme (San Francisco: ASP), 19 

Kotani, T., Kawai, N., Matsuoka, M. and Brinkmann, W. 1996, PASJ, 48, 619 

Kubota K., et al. 2010, ApJ, 709, 1374 

Laurent, P. k Titarchuk, L. 1999, ApJ, 511, 289 

Laurent, P. k Titarchuk, L. 2010, submitted to ApJ 

Laurent, P. k Titarchuk, L. 2007, ApJ, 656, 1056 

bibitem [Laurent k Titarchuk (1999)]LT99 Laurent, P. k Titarchuk, L. 1999, ApJ, 
511, 289 

Lyubarskii, Yu. 1997, MNRAS, 292, 679 
Margon, B. 1984, ARAA, 22, 5072 

Marshall, H. L., Canizarcs, C. R., Schulz, N. S., 2002, ApJ 564, 941 
Miller, J, 2007, AR&A, 45, 441 

Munoz-Darias, T., Casares, J., k Martinez-Pais, I. G. 2008, MNRAS, 385, 2205 
Nandi, A., Chakrabarti, S.. Belloni, T., Goldoni, P. 2005, MNRAS. 359, 629 
Revuivtsev, M. et al, 2006, AVA, 147. 545 



20 


Romney, I. D., Schilizzi, R. T., Fejet, I., Spencer, R. E., 1987, ApJ, 321, 822 

Safi-Harb, S., k Kotani, T. 2003, in New Views on Microqnasars, ed. P. Durouchoux, Y. 
Fnchs, & J. Rodriguez (Kolka.ta: Centre for Space Physics), 279 

Seifina, E. 2000, Astron. Rep., 43, 305 

Shaposhnikov, N., k Titarcliuk, L. 2009, ApJ, 699. 453 (ST09) 

Sunyaev, R.A. V Titarchuk, L. 1980, AkA, 86, 121 
Titarchuk, L. & Seifina, E. 2009, ApJ, 706, 1463 (TS09) 

Titarchuk, L. , Shaposhnikov, N. k Aref’ev, V. 2007, ApJ, 660, 556 
Titarchuk, L. k Shaposhnikov, N. 2005, ApJ, 626, 298 
Titarchuk, L. k Zannias, T. 1998, ApJ, 499, 315 
Titarchuk, L. , Lapidus, 1. 1. k Muslimov, A. 1998, ApJ, 499, 315 
Titarchuk, L. , Mastichiadis, A k Kylafis, N. 1997, ApJ, 487, 834 

van der Klis, M. 1995, in X-ray binaries, ed. W.H.G. Lewin, J.,van Paradijs, k E.P.J. van 
der Heuvel (Cambridge Univ. Press, Cambridge), p. 252 

Zwitter, T., Calvani, M., k D’Odorico, S., 1991, AkA, 251, 92 


This preprint was prepared with the AAS PTqX macros v5.2. 



Table 1. 


The list of sets (groups) of RXTE observation of SS 433 


Number of set 

Dates, MJD 

RXTE Proposal ID 

UT Dates 

Type of Light Curve 

Ref. 

R1 

50191-50194 

10127 

Apr. 18 - 21, 1996 


1 

R2 

50868-50907 

20102, 30273 

Feb. 24 - Apr. 4, 1998 

outburst 

1,2 

R3 

52222-52238 

60058 

Nov. 9 - 11, 2001 

outburst 

1 

R4 

52544.46-52544.74, 52913-52914 

70416, 80429 

Apr. 18 - 21. 2002 


1 

R.5 

53076-53092 

90401 

Mar. 12 - 28, 2004 

outburst 

1, 2 

R6 

53239-53610 

90401, 91103, 91092 

Jul. 28 - Aug. 28, 2005 

outburst decay 

this work 

RT 

54085-54096 

92424 

Dec. 17 - 27, 2006 

outburst decay 

this work 


References: (1) Filippova et al. (2006), (2) Nandi et al. (2005). 


Table 2. The list of GRS 1915+105 observations used in analysis. 


Satellite 

Obs. ID 

Start time (UT) 

End time (UT) 

BeppoSAX 

20985001 1 

2000 Apr. 21 08:55:30 

2000 Apr. 21 15:16:47 



Table 3. Best-fit parameters of spectral analysis of PCA and HEXTE observation of SS 433 in 3-150 keV energy 

range*. Parameter errors correspond to la confidence level. 
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60058-01-09-00* 52230.30 0.717 0.7576 1.30(5) 0.5(2 ) 2.1(1) 6.92(4) 0.73(6) 2.18(3) 8.69(4) 0.25(8) 0.79 (69) 3.46/2.17 

60058-01-10-00* 52232.15 0.859 0.7690 1.2(1) 2.00 2.27(3) 6.96(3) 0.86(4) 3.2(3) - - 1.2 (69) 4.15/2.74 

60058-01-11-00* 52233.27 0.952 0.7758 1.20(8) 2.00 2.06(1) 6.97(4) 0.63(7) 1.9(2) 8.7(2) 0.2(1) 1.2 (69) 3.57/2.37 
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91103-01-02-01 53580.83 0.952 0.0865 1.18(3) 2.00 1.88(5) 7.1(1) 0.78(2) 1.7(5) 8.7(1) 0.2(2) 1.18(73) 3.38/2.13 

91103-01-02-00* 53580.90 0.958 0.0870 1.18(3) 0.2(2) 2.21(2) 6.70(2) 0.27(3) 0.6(3) 8.0(2) 0.4(1) 0.97(69) 2.97/1.74 

91103-01-03-00 53581.53 0.006 0.0908 0.94(6) 0.58(5) 1.6(1) 6.96(3) 0.99(3) 1.61(6) 9.3(1) 0.10(8) 1.02(73) 2.79/2.30 
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91092 - 02 - 07-00 53597.34 0.214 0.1883 1 . 20 ( 3 ) 2.00 3 . 1 ( 2 ) 6 . 88 ( 2 ) 0 . 46 ( 4 ) 2 . 2 ( 2 ) 8 . 5 ( 1 ) 0 . 17 ( 8 ) 1.10 ( 72 ) 5 . 41 / 3.73 

91092 - 02 - 08 - 00 * 53598.38 0.293 0.1947 1 . 20 ( 1 ) 0 . 1 ( 1 ) 2 . 64 ( 9 ) 6 . 80 ( 6 ) 0 . 53 ( 2 ) 1 . 6 ( 3 ) 8 . 3 ( 2 ) 0 . 2 ( 1 ) 1.13 ( 71 ) 3 . 358 / 1.71 

91103 - 01 - 10-00 53610.77 0.241 0.2712 1 . 18 ( 4 ) 2.00 2 . 59 ( 3 ) 6 . 63 ( 2 ) 0 . 40 ( 4 ) 2 . 3 ( 1 ) - - 1.08 ( 73 ) 4 . 12 / 2.75 
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when parameter log(A) > 1, then it is fixed at 2.0 (see comments in the text), tt normalization in units of 10 3 ( 10 10 ergs/s/cm 2 ), spectral flux in units of 

10~ 10 erg/s/crn 2 , au ne 2 for Gaussian2 component was fixed at 0.01 keV. 
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Table 5. Best-fit parameters of spectral analysis of PCA and HEXTE observation of 
SS 433 in 3-150 keV energy range in the model: wabs*( bmc+gaussian+gaussian+ ,, bbody”) 
for observations with numbers 91103-NN-NN-NN. Parameter errors correspond to ler 

confidence level. 


Model 

Parameter 

-01-02-00 

-01-05-00 

-01-05-01 

-01-06-00 

-01-06-01 

-01-07-00 

-01-08-00 

-01-09-01 

bmc 

r 

2.18(3) 

2.2(1) 

2.20(6) 

2.19(5) 

2.2(1) 

2.21(1) 

2.2(1) 

2.2(1) 


kT (keV) 

1.04(1) 

1.29(8) 

1.25(6) 

1.22(7) 

1.23(6) 

1.21(4) 

1.18(7) 

1.11(9) 


logA^ 

0.2(2) 

0.01(4) 

0.9(1) 

2.0 

2.0 

2.0 

2.0 

2.0 


KL 

2.21(2) 

2-5(3) 

3.1(4) 

3.2(3) 

3.2(3) 

3.3(2) 

3.0(4) 

3.0(4) 

Gaussian 1 

Efinel (keV) 

6.70(2) 

7.08(6) 

7.02(7) 

7.08(5) 

7.06(3) 

7.03(3) 

7.12(9) 

6.87(6) 


mine 1 (keV) 

0.27(3) 

0.75(8) 

0.75(9) 

0.90(9) 

0.75(5) 

0.75(6) 

0.98(6) 

0.84(9) 


Np 

' Uriel 

0.6(3) 

2.8(5) 

3.0(5) 

3.0(4) 

3.5(3) 

2.8(3) 

4.1(3) 

3.9(7) 

Gaussian2 

E line.2 (keV) 

8.0(2) 

_ 


7.0(1) 

7.01(8) 

7.1(4) 

7.0(1) 

. 


wtt 
iN line 2 

0.4(1) 

- 


0.05(3) 

0.6(3) 

0.9(2) 

0.5(2) 

- 

’'bbody 1 ' 

f "bbody" 

(keV) 

4.9(2) 

4.9(2) 

4.9(4) 

5.1(9) 

5.1(3) 

5.6(5) 

5.2(4) 

5.5(3) 


N +t 

bbody” 

1.17(7) 

2.1(6) 

1.2(8) 

0.8(6) 

2.6(5) 

1.7(4) 

1.9(4) 

2.3(6) 


EW” bbody ” 

1.2(2) 

1.3(3) 

1.4(5) 

1.4(5) 

1.4(1) 

1.4(2) 

1.4(1) 

1.4(2) 


(koV) 









Fluxttt 

3 - 60 keV 

4.99 

3.90 

5.41 

5.56 

6.72 

6.30 

6.56 

6.54 


13 - 150 keV 

3.17 

2.13 

3.34 

3.44 

4.01 

3.88 

4.29 

4.34 


X 2 (d.o.f.) 

1.17(69) 

0.97 (69) 

1.05 (69) 

1.25 (69) 

1.01 (69) 

1.04(69) 

1.03 (69) 

1.05(69) 


f when parameter log(A) > 1, it is fixed at 2.0 (see comments in the text), ^ normalization in units of 10 '*(10” 10 
ergs/s/cm 2 ), ttt spectral flux in units of 10“ 10 erg/s/cm 2 , <7j ine 2 for Gaussian2 component was fixed at 0.01 keV. 
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Table 6. Best-fit parameters of spectral analysis of PCA and HEXTE observation of 
SS 433 in 3-150 keV energy range in the model: wabs*( bmc+gaussian+gaussian+”bbody”) 
for observations with numbers 91092-NN-NN-NN. Parameter errors correspond to la 

confidence level. 


Model 

Parameter 

-01-01-00 

-01-02-00 

-02-04-00G 

-02-05-00 

-02-06-00 

-02-06-01 

-02-08-00 

bmc 

r 

2.20(3) 

2.20(3) 

2.22(3) 

2.2(3) 

2.22(1) 

2.22(3) 

2.20(1) 


kT (keV) 

1.1(1) 

1.22(3) 

1.20(6) 

1.28(6) 

1.21(7) 

1.3(1) 

1.35(2) 


logA* 

0.6(3) 

0.1(1) 

0.6(4) 

0-5(3) 

1-2(1) 

0.6(2) 

0.1(1) 


Np 

bmc 

2.71(2) 

2.90(5) 

2.7(1) 

2.9(1) 

2.6(2) 

2.67(9) 

2.64(9) 

Gaussian 1 

E,inel (keV) 

6.89(3) 

6.88(1) 

6.81(2) 

6.83(1) 

6.90(2) 

6.90(2) 

6.80(6) 


ariinel (keV) 

0.75(6) 

0.50(2) 

0.48(5) 

0.46(4) 

0.51(5) 

0.53(3) 

0.53(2) 


M + t 
'line 1 

3.7(4) 

2.8(1) 

1.7(1) 

2.6(2) 

2.4(2) 

2.5(1) 

1.6(3) 

Gaussian2 

E, ine2 (keV) 

9.1(1) 

. 

8.37(5) 

8.44(1) 

8.4(9) 

8.50(8) 

8.3(2) 


Mtt 

Une2 

0.2(1) 

- 

0.40(7) 

1.2(4) 

2.3(9) 

0.32(7) 

0.21(2) 

"bbody” 

F” bbody" 

4.5(2) 

4.5(2) 

4.4(4) 

4.9(6) 

5.5(5) 

5.0(4) 

4.5(4) 


(keV) 









N ” 66or/i,” 

1.7(7) 

2.2(4) 

1.3(4) 

2.8(5) 

1.9(6) 

1.7(1) 

1.9(3) 


EW M 5fj 0{ fy” 

1.2(2) 

1.3(3) 

1-4(5) 

1.4(5) 

1.4(1) 

1.4(2) 

1.4(1) 


(keV) 








Fluxttt 

3-60 keV 

4.36 

4.07 

4.42 

4.43 

3.61 

3.98 

3.35 


13 - 150 keV 

2.56 

2.44 

2.80 

2.90 

1.46 

1.81 

1.71 


X 2 (d.o.f.) 

0.96(69) 

1.07 (69) 

1.39 (69) 

1.06 (69) 

1.05 (69) 

0.96(69) 

1.05 (69) 


t when parameter log(A) > 1, it is fixed at 2.0 (see comments in the text), tt normalization in units of 10 3 (10 10 
ergs/s/cm 2 ) x 10"' 3 , ttt spectral flux in units of 10 _to erg/s/cm 2 , n Une2 for Gaussian2 component is fixed at 0.01 keV. 


Table 7. BH masses and distances 


Source 

Kiyn 

(M©) 

i orb 

(deg) 

l scal 

(deg) 

D h 

(kpc) 

M.sca/ 

(M 0 ) 

fiscal 

(kpc) 

Refs 

GX 339-4 

>6 


70 

7.5±1.6 

12.3±1.4 

5.75±0.8 

1. 2 

SS 433 

4.3±0.S 

78.7 

80 

5.0±0.5 

2 2 


3, 4, 5 


“ Dynamically determined BH mass, b Source distance found in literature. 

References: (1) Munoz-Dias et al. (2008), (2) Hynes et al. (2004), (3) Romney et al. (1987), (4) Hillwig k. Gies (2008), (5) 

M argon (1984). 




Converging inflow Bulk comptonization 


Fig. 1. Schematic view of the proposed geometry for thermal and bulk Comptonization 
regions in the source SS 433 hosting a BH with PL-like emission at high energies. The 
terminal direction of viewing of the disk central part to Earth observer marks by vector l. 
Here vector n is denoted normal to orbital plane, i is an inclination angle, ip is a precession 
phase angle, a is cone precession angle and a bended arrow points the direction of the disc 
precession motion around normal n to the orbital plane. 
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Fig. 2.-— Two representative spectra of SS 433 in EF(E) units for outburst IS (top) and for 
quiet IS (bottom). Blue line presents XSPEC continuum component of the best-lit model (sec 
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Fig. 3. The best-fit spectrum during of 2004 outburst transition in EF(E) units {top) 
and in normalized counts units {bottom, panel ) with x for the 90401-01-01-01 observation 
(R5 set). Left: the best fit which does not include a high-temperature “bbody” component 
(Xred = 4-05 for 72 d.o.f.) and right: the best-fit spectrum and Ay 2 , when residual hump at 
20 keV is modelled by a high-temperature “bbody” component with Tbb = 5.9 ± 0.2 keV 
(xiw—l-OO for 70 d.o.f.). On top panel data are denoted by black points. Model consists of 
four components shown by blue, dashed purple, light blue, and red lines for BMC, Gaussianl , 
GaussianS and bbody components respectively. 
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Fig. 4.-— The best-fit spectrum during of 2005 outburst transition in EF(E ) units (top) and 
in normalized counts units ( bottom panel) with \ for the 91092-01-01-00 observation (R6). 
Left : the best fit which does not include a high-temperature “bbody” component (M? efl >=5.01 
for 72 d.o.f.) and right: the best-fit spectrum and Ay 2 , when residual hamp around 14 keV is 
modelled by a high-temperature bbody” component with Tbb — 4.5 ±0.2 keV (Xred = 1-07 for 
70 d.o.f.). On top panel data are denoted by black points. Model consists of four components 
shown by blue, dashed purple, light blue, and red lines for BMC , Gaussianl, GaussianS and 
bbody components respectively. 
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Fig. 5.— From top to bottom: Evolution of the; flux density (GBI, 2.25 GHz - blue, 8.3 GHz 
- green), RXTE/ASM count rate, spectral index a ra dio in corresponding radio GBI range, 
BMC normalization and photon index F for the middle of 2004 outburst of SS 433. Blue 
strips mark eclipsed intervals of light curves (around MJD 50872-50875, 50886-50889 and 
50899-50 
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Fig. 6. — From top to bottom: Evolution of the flux density [0.96 (green) and 2.3 (blue) 
Gz), RATAN-600], RXTE/ASM count rate, BMC normalization, photon index F and color 
temperature Tbb of high temperature “Black Body” spectral component for the middle of 
2001 outburst of SS 433. Blue strips mark eclipsed intervals of light curve around MJD 
52219 52221 and 522222 52234 during R3 set. 
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Fig. 7.— From top to bottom: Evolution of the flux density (2.3 GHz, RATAN-600), 
RXTE/ASM count rate, BMC normalization and photon index F for the middle of 2004 
outburst of SS 433. Blue strips mark eclipsed intervals of light curve around MJD 53070 - 
53100 (R5 set). 
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Fig. 8. — Photon index T plotted versus BMC normalization (left) and Comptonized frac- 
tion (right) for all observations. Blue squars correspond to 2004 outburst. Crimson circle 
marks 2005- 2006 observations. Red squares is selected for decay of 2005 events. The rest 
observations denote by black triangles. 
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Fig. 9. — From top to bottom: Photon index T versus orbital phase (p during MJD=50191- 
50194 (El), 50897-50907 (E2), 52225-52238 (E3), 53579-53588 (E4). Blue vertical strip 
marks an interval of the primary (optical) eclipse (Goranskij et al. 1998). Points marked 
with rose oreol correspond to spectra fitted by the model which includes ’’high- temperature 
BB” component (see also Fig. 1 and Tables 4-6). Red points ( for two lower panels) correspond 
to observations during the primary eclipse when the ” high-temperature BB” component is 
not detected, although it presents in spectral residual before and after this eclipse. 
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Fig. 11. Top : From up to bottom: Optical V-band light curve, BMC normalization and 
photon index F for RG set. Points A and B mark moments 53581 and 53585 (at eclipse 
and after eclipse) respectively. Blue strips mark eclipsed intervals of light curve. Bottom-. 
PDSs ( left panel ) are plotted along with energy spectral diagram E x F(E) (right panel ) for 
observational points A ( 91103-01-03-00 , blue) and B (91103-01-07-00, red). 
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Fig. 14. Broad-band BeppoSAX spectrum of GRS 1915+105 (top) and RXTE/PCA and 
IBIS/ISGRI/INTEGRAL spectrum of SS 433 ( bottom ) in EF(E) units obtained during in- 
termediate transition state of GRS 1915+105 on April 21, 2000 (ID=209850011) and simul- 
taneous RXTEj INTEGRAL observations of SS 433 for outburst transition state on March 
24 27. 2004. 




